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Abstract 

We survey the close relationship between sound and light rays and geometry. In the case 
where the medium is at rest, the geometry is the classical geometry of Riemann. In the 
case where the medium is moving, the more general geometry known as Finsler geometry is 
needed. We develop these geometries ab initio, with examples, and in particular show how 
sound rays in a stratified atmosphere with a wind can be mapped to a problem of circles 
and straight lines. 



1 Introduction: Propaga- 
tion of Waves in Moving 
Media 

Almost everyone must be familiar, by report at 
least, with optical mirages. In hot regions such 
as the desert, when the speed of light is greater 
at low altitudes than at higher altitudes, dis- 
tant palm trees can appear inverted as if re- 
flected off a cool pool of water in a nearby oasis, 
see Figure [I] It is less well known that in cold 
regions such as the arctic, where the opposite 
conditions prevail, distant ships can appear in- 
verted in the sky, the light from them having 
been bent over an intervening iceberg, see Fig- 
ure [2j More mundanely, motorway travellers 
on hot dry summer days often have the discon- 
certing impression that there are sheets of wa- 
ter lying on the road some distance ahead. The 
explanation of phenomena like this is easily un- 
derstood using the concept of light rays subject 
to Snell's Law for a stratified medium. Alter- 
natively we may apply Huygens's wave theory 



according to which, in passing over the iceberg, 
the higher parts of the wave front move faster 
than the lower parts causing it to bend down- 
wards. 

Mutatis Mutandis similar phenomenon can 
apply to sound waves whose speed increases 
with increasing absolute temperature above 
zero, T, as J 1 '. During a warm sunny day 
therefore, when the temperature is typically 
hotter near the ground due to the sun's rays, 
the sound waves should be bent upwards. 
On a clear night, when the temperature near 
the ground drops rapidly by radiative cooling, 
sound waves should be bent downwards allow- 
ing distant sources, such as cars on a motorway, 
to be heard much more clearly than during the 
day. One of us, living as he does a kilometre or 
so away from the A14, which is, judged by the 
amount of traffic it carries effectively a motor- 
way, had until recently always supposed that it 
was this temperature effect that was responsi- 
ble for the din experienced at times when con- 
templating the night sky in his garden. How- 
ever, the temperature effect should be isotropic, 
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Figure 1: A 19 th century woodcut showing 
a mirage in the desert. Source: 'Elements de 
Physique ', Ganot 



in other words it should affect the noise of cars 
coming from all directions and including those 
on less busy but closer city roads. The greater 
volume and speed of traffic on the A14 alone 
should not be so overwhelming. 

The obvious directional influence is that of 
the wind. However, the speed of sound v s rs 
1250 km per hour is so much greater than typi- 
cal wind speeds v « 30 km per hour that simple 
convection of sound waves by the wind cannot 
be responsible for any significant directional ef- 
fect. As pointed out to one of us by a colleague, 
Hugh Hunt, who made a similar point in the 
New Scientist of 15th April 2009 in response to 
readers' queries, it is not the wind velocity, but 
its gradient, that is its variation with height, 
called wind shear or more technically vorticity 
which is important. As is clear from observing 
clouds, wind speeds, while remaining roughly 
horizontal and much slower than the speed of 
sound, increase quite sharply with height, while 
remaining roughly horizontal. Thus it is not 
just the increased velocity |v| that matters, but 
its gradient or vorticity uj — curl v. 

Of course, examination of the literature 
shows that this is not a new observation and 
many papers and textbooks contain a simple 



Figure 2: Etching of a sketch by the explorer 
William Scoresby, showing an arctic mirage. 
Source: 'The life of William Scoresby', Scoresby- 
Jackson 

qualitive discussion. The earliest we have dis- 
covered dates from 1857 [T| and is due George 
Gabriel Stokes (1819-1903) who was appointed 
to the Lucasian Chair of Mathematics in Cam- 
bridge in 1849 and held it until his death 54 
years later. His explanation, elaborated on by 
Osborne Reynolds (1842-1912) in 1874 0, fol- 
lowed very closely Huygens's explanation for re- 
fraction by a gradient in the refractive index. If 
the wind direction is towards us, and the wind 
speed increases with height, the higher parts of 
the wave front move faster than the lower parts 
and the wave front is bent over towards us. If 
on the other hand the wind direction is in the 
opposite direction, the wave front will be bent 
upwards and hence away from us. The net ef- 
fect is that if the vorticity u> is non zero, the 
sound rays are deflected in an analogous way 
to the deflection of a charged particle of mass 
m and charge e by the Lorentz Force due to 
magnetic field B. In fact this is not just an 
analogy but, as we shall see shortly, a precise 
correspondence for low wind speeds: 

w = -— B. (1) 

m 

Note that since it is the velocity gradient that 
matters, it is not the direction of the wind at 
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Figure 3: The cycloid, as traced out by a point on the circumference of a rolling circle 



the ground or at a height above , but their dif- 
ference which is important. In practice how- 
ever, the speed of the wind is almost always 
much slower near the ground than above, and 
so usually it is the velocity at higher altitudes 
which matters. Thus, in principle there are two 
effects acting, and which is the more important 
depends on which induces the greater curvature 
K to the sound rays. As Stokes and Reynolds 
realised n ~ i|r = hh'rf- for the thermal ef- 

n oz 2 T Oz 

feet and k ~ — for the wind shear, where 

v s oz 1 

n is the index of refraction. It is probably no 
coincidence that Stokes's work followed shortly 
after the first demonstration, by the German 
physicist Karl Friedrich Julius Sondhauss (1815 
- 1886) in 1853 [3], that a balloon filled with 
CO2 will act as a sound lens, focussing the 
sound of a ticking watch so as to render it au- 
dible some distance away. 

Examples of the interplay of these two ef- 
fects, which can cause sound to travel over large 
distances, abound. In early June, 1666, during 
the war between the Dutch and the English, 
both Samuel Pepys and John Evelyn reported 
in their diaries that while the sound of gun fire 
from ships off the coast of Kent could be heard 
clearly in London, it was not audible at all at 
the ports of Deal or Dover. As Pepys at the 
time observed: 

"This . . . makes room for a great 
dispute in Philosophy: how we 
should hear it and not they, the 
same wind that brought it to us be- 
ing the same that should bring it to 
them." 

Infrasound (sound of very low frequency) from 
the volcanic explosion of Krakatoa on Au- 
gust 27 1883 was heard to travel several times 
around the earth. During the first world war, 



the noise of the very large guns on the Western 
Front was often audible within a range of a 100 
km or so, and often beyond 200 km, but not 
within a "zone of silence" between 100 and 200 
km. On September 21, 1921 there occured an 
enormous explosion at Oppau on the Rhine and 
the same phenomenon was observed, see Figure 
® 

One explanation for some of these phenom- 
ena is reflection of sound from a layer of air in 
the upper atmosphere with a higher sound ve- 
locity. The dependence of the velocity of sound 
in the atmosphere follows its temperature pro- 
file. In the absence of wind, the temperature 
(and hence velocity of sound) decreases under 
normal circumstances up to a height of about 
10 km, remains constant up to roughly 25 km, 
then incresases to around 50 km, where it has 
a local maximum, and then has a local mini- 
mum at about 80 km, see figure [5j A simple 
application of the law of refraction 

n{z) sm9(z) = constant, (2) 

where the local speed of sound v s (z) cx ^jy, re- 
veals that sound rays can bounce off the max- 
ima of v s (z) but can also be trapped at the 
minima. Presumably the latter effect was re- 
sponsible for the long distance propagation of 
infrasound from Krakatoa. A similar phenom- 
ena occurs for sound waves in the ocean. The 
speed of sound initially decreases with depth 
and then increases exhibiting a minimum value 
at a depth of around 1 km at the so-called SO- 
FAR channel. This allows whales to communi- 
cate over very large distances [8j |9] and more 
sinisterly, submarines to snoop on other sub- 
marines. 

For simple profiles n{z), the law of refrac- 
tion ([2| allows simple solutions for the rays. 
Thus if n{z) cx z, one has catenaries with the 
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Figure 4: The explosion at Oppau, 1921. 
Filled circles represent locations where the ex- 
plosion was heard and empty circles locations 
where nothing was heard 

Source: 'Strange sounds in the Atmosphere: Part I', 



Figure 5: The variation of temperature and 
sound speed with altitude 



R. K. Cook, Sound 1, 12-16 (1962). 



horizontal axis z — as the directrix. This 
gives a rough description of mirages in the 
desert. If n(z) oc ^ , the rays are cycloids, the 
curve traced out by the point on the circumfer- 
ence of a circle rolling without slipping along 
the horizontal axis (see Figure |3j. One can 
imagine this as the path of a glow-worm sitting 
on the rim of a bicycle wheel as it rolls along in 
the dark. This could decsribe the rays passing 
over icebergs in the arctic. More importantly 
for what follows later, this can also be achieved 
by assuming that n(z) oc -, in which case the 
rays are semi-circles centred on the horizontal 
axis. The addition of the effects of wind com- 
plicates considerably this simple picture. 

Considering for a moment fundamental 
physics, one of the clearest trends in research 
over the last 100 years or so has been what one 
might call the geometrisation of physics. This 
is by no means a new phenomenon. Plato's as- 
sociation of the regular, or Platonic, solids with 
the classical elements is perhaps the earliest at- 

1 Throughout this paper we will use the convention 
summed over 



tempt to explain the world through geometrical 
intuition. Geometry, in the modern sense of 
diffcrcntiable manifolds, really caught hold in 
physics in 1915 with Einstein's general theory 
of relativity. Since then, the use of geometry 
has accelerated. Modern string theory, with 
its 10 space-time dimensions and complicated 
internal 6-dimensional Calabi-Yau manifolds is 
perhaps the clearest example of this geometri- 
sation process. 

A geometrical approach has much to recom- 
mend itself, even to describe the more concrete 
physics we have thus far been considering. In 
this article, we will discuss how the properties 
of sound and light rays can be considered in a 
geometrical light. We will be interested in the 
behaviour of solutions of the wave equation^] in 
a moving medium: 

'(?r- W 1 ^-) 2 - C 2 h ij = 9 l 1 U (x, t)=Q. 

Adt dx l ) dx l dxJ\ v ' ; 

(3) 

hat indices i, j = 1 . . . n which are repeated should be 
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Here h %J is related to the speed of sound in 
a direction parallel to the unit vector by 
v m = c(h % ' j ' minij)^ 1 / 2 and W l is a vector giving 
the velocity of the medium at each point. We 
will be interested in particular in disturbances 
with short wavelength, which move along rays. 
We will start by considering the case of a static 
medium, W l — 0. The geometry of the rays in 
this case is well known and we shall describe 
some of its properties and give some examples. 
We will then move on to the case of a mov- 
ing medium and demonstrate that this can also 
be described geometrically, as we showed in [5] . 
We will have to loosen slightly our notion of 
a geometry in order to do so but the type of 
geometry which arises, a Finsler geometry, is 
very natural. In fact, this geometry first arose 
in our studies of light rays near rotating black 
holes. The link between the effects of gravity 
on light rays and refraction of light and sound 
waves by media can be made explicit and is 
the basis for much current work on optical and 
acoustic black holes. These analogue models 
for black holes allow experimentation in a lab- 
oratory, which would of course not be possible 
with real black holes. For this reason, a thor- 
ough understanding not just of the mechanisms 
of refraction, but the geometry of refraction is 
of great relevance both for terrestrial and celes- 
tial physics. 

2 Geodesies and Fermat's 
Principle 

2.1 Optical and Acoustic Geome- 
tries 

The elementary theory of mirrors and lenses 
is largely concerned with tracing the paths of 
light rays on reflection (Catoptrics) or refrac- 
tion (Dioptrics) at a surface. Although the 
Greeks were uncertain whether light proceeds 
from the eye to the object (emission theory) or 
from the object to the eye intromission theory, 
and whether its passage is instantaneous or at a 
finite speed, nevertheless Heron of Alexandria 
(c. 10-70 AD) was able to formulate the laws 
of reflection in terms of a Principle of Short- 



est Length from object to eye via the reflecting 
surface. The occurrence of caustics shows that 
there can be more than one path, not all of 
which are necessarily the shortest and more ac- 
curately we refer to the principle Principle of 
Stationary Length, that is the length of each 
ray is merely stationary among all neighbour- 
ing paths. 

Despite important pioneering work by Ibn 
al-Haytham or Alhazen (965-C.1040) demolish- 
ing the emission theory and investigating re- 
fraction, it took longer to unravel the funda- 
mental law of Dioptrics. It it was not un- 
til independent work by Abu Sa'd al-'Ala' ibn 
Sahl (c.940-c.l000), Thomas Harriot (c.1560- 
1621) and Willebrord Snel Van Royen (1580- 
1626) that the familiar Law of Sines was es- 
tablished and Pierre Fermat (c 1601- 1665) for- 
mulated his unified Fermat's Principle of Sta- 
tionary Time. The idea is that the slowness of 
the ray inside a medium is proportional to its 
refractive index n. Note that the finite speed 
of light was only demonstrated by by Ole Roe- 
mer (1644-1710) using the eclipse of Jupiter's 
moon Io in 1676. The relation of the refractive 
index to the speed of light remained contro- 
versial until experiments by Armand Hippolytc 
Louis Fizeau (1819-1896), Fresnel Augustin- 
Jean Fresnel (1788-1827) and George Biddell 
Airy (1801-1892) in the nineteenth century fi- 
nally established its velocity as ~, where c is its 
speed in vacuo. According to the discredited 
corpuscular theory the opposite relation holds. 
Both theories give the same rays but the speed 
with which light follows the rays differs. It 
would be more accurate therefore to speak of 
Fermat's Principle of Stationary Optical Path 
Length , where the optical length of a path 7 is 
given by 

L= [ n(x)|dx|, (4) 

where we have allowed for the possibility that 
the refractive index may depend upon position, 
as for example it docs in a vertical stratifed 
medium such as we encounter discussing mi- 
rages. Fermat's principle becomes the Varia- 
tional Principle 

5L = 5 ( n(x)\dx\ =0. (5) 
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By the time Fermat introduced this principle, 
Christiaan Huygens (1629-1695) had initiated 
his wave theory of light and derived Fermat 's 
Principle from it. His derivation makes it clear 
that it is the optical length or optical distance 
which enters into all interference effects, and 
its therefore appropriate to say that all opti- 
cal measurements measure Optical Geometry. 
By the same token, measurements using sound 
waves may be said to measure Acoustic Geom- 
etry and measurements using seismic waves, as 
on the earth or more recently the moon [TU] to 
measure Seismic Geometry. 

It is clear that these geometries will, in an 
inhomogcneous medium for which the depen- 
dence of n on x is non-trivial, differ consid- 
erably from Euclidean Geometry. The exis- 
tence of such non-Euclidean Geometries was 
first realised by pure mathematicians in the 
early part of the nineteenth century working 
on the foundations of geometry. For centuries, 
people had been attempting to derive, start- 
ing from the other aximons, Euclid's fifth ax- 
iom: that through any point not on a given line 
there is exactly one line parallel to the first. 
This seemed to them so obvious that it was 
"neccessarily" true. Eventually they gave up, 
Johann Carl Friedrich Gauss (1777-1855) pri- 
vately and Janos Bolyai (1802 -1860) and Niko- 
lai Ivanovich Lobachevsky (1792-1856) pub- 
licly, showed the existence of two other types 
of homogeneous and isotropic Congruence Ge- 
ometries, Spherical and Hyperbolic. The first 
is easy to grasp in two dimensions since it is 
just the geometry of the standard sphere, S 2 
with the stationary paths (or geodesies) being 
the great circles. Navigators, either by sea or 
air, have been using spherical geometry since 
at least the time of Columbus. It is not too 
difficult to imagine a sphere in one dimension 
higher and indeed if the refractive index were 
to vary as 

n — t. — 

the optical metric would be precisely that of 
three dimensional spherical space S 3 of radius 
UqR. The resulting optical device is known as 
Maxwell's Fish Eye Lens since all rays emanat- 
ing from any point x e are circles which recon- 



verge onto the antipodal point x e = t^tj- Wc 
can verify that the optical distance along a ra- 
dial geodesic is given by: 



(G) 



n 



1 



-dr = nnoR, 



(7) 



4R 2 



which is finite. Because, for large enough |x|, 
the refractive index n drops below unity, the 
construction of such a device would require the 
manufacture of a suitable "meta-material" . A 
more practical device was invented by Rudolph 
Karl Luneburg (1903-1949) and has 

n = (2|x| 2 -l)), |x|<l; 

n = 1, |x|>l. (8) 

This will focus all rays incident on it a fixed 
direction to the point on the circumference in 
the opposite direction. 

The radius of curvature of the Spherical 
Geometry given by (p| is n a R. To obtain 
Hyperbolic Space H , often called Bolyai- 
Lobachevsky space or just Lobachevsky Space, 
one needs only to let the radius of curvature 
become pure imaginary which leads to 

l 1 AR 2 > 

The refractive index becomes infinite when 
|x| = 2R which should be thought of as the 
boundary at 'infinity' of hyperbolic space. In 
fact the reader may easily verify that the opti- 
cal distance along a radial path from the origin 
x = to the boundary at infinity is, by contrast 
with the case of spherical space, infinite. 

Jules Henri Poincare (1854-1912) gave a 
simple analogy which illuminates the roles 
played by the flat Euclidean metric geometry 
for which n = 1 and the curved non-Euclidean 
geometry for which n is given by ([9} as fol- 
lows. Imagine a medium whose temperature 
varies with radius as — occupying a ball of ra- 
dius |x| = 2R, as measured by a measuring 
rod made from a substance such as Invar whose 
length is independent of temperature. If mea- 
sured by a different ruler made of a material 
which expands in proportion to the tempera- 
ture, then the rod will shrink to zero length as 
it approaches the the boundary which is at zero 
temperature. The boundary will thus seem to 
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be infinitely far away as measured by the sec- 
ond measuring rod. At the time that Poincare 
wrote, people were still reeling under the dis- 
covery that what seemed to have been well es- 
tablished since antiquity: that Euclid's Geo- 
metrical Axioms were not logical necessities. 
There was thus great interest among geometers 
and philosophers as what was the "correct" or 
"real" geometry of space. For Poincare it all 
depended upon how you measure it. He did 
however believe that one should always be able 
to find a system of measurements in which Eu- 
clid's Geometry holds. 

Despite appearances both Spherical Geom- 
etry and Hyperbolic Geometry are, like their 
simpler Euclidean counterparts, both isotropic 
and homogeneous. For this reason they are can- 
didates for the physical geometry of space, as 
measured for example by light rays in vacuo. 
Indeed according to Einstein's theory of Gen- 
eral Relativity just these three possibilities 
can arise in the theory of the Expanding Uni- 
verse proposed by Alexander Alexandrovich 
Friedman (1888 -1925) and Monsignor Georges 
Henri Joseph Edouard Lemaitre (1894 - 1966). 
For many years cosmologists have been at- 
tempting to decide which best fits the observed 
universe. Based on observations of the Cos- 
mic Microwave Background (CMB), and other 
data, the consensus now is that it is flat Eu- 
clidean geometry. 

Of course no realistic medium is exactly ho- 
mogeneous or isotropic, and in particular the 
speed of light, or sound, may depend upon di- 
rection as well as position. A familiar example 
is provided by bi-refringence in crystals such 
as calcite for which the ordinary and extra or- 
dinary ray have different refractive indices, n a 
and n e . This more general situation may be 
taken into account using a more general geom- 
etry invented by Georg Friedrich Bernhard Ric- 
mann (1826 - 1866) called Riemannian Geome- 
try in which the optical or acoustic path is given 

by 



SL 



dx l dxi 



0, (10) 



where hij is, in three spatial dimensions, a 3 x 3 
symmetric array called the metric tensor. In- 



deed one frequently introduces what is called a 
line element which is an expression for the in- 
finitesimal form of a generalised Pythagoras's 
theorem: the infinitesimal distance ds between 
points x l and x l + dx l is given by 



ds 2 



hij{xk)dx l dx 3 



(11) 



Thus, for example, in the case of a uniaxial bi- 
refringent medium with unit field n in the ordi- 
nary direction the Joets-Ribotta optical metric 

is a 



ds 2 = n^dx 2 + (n 2 , - rig)(n.(ix) 2 



(12) 



Ricmann's ideas are fully incorporated into 
Einstein's Theory of General Relativity. In- 
deed, for a static spacetime a simple form of 
Fermat's Principle holds which, for example, 
allows one to discuss the optics of black holes 
in terms of an effective refractive index (in so- 
called isotropic coordinates) 



n(x) 



GM 



2c 2 



1 \° / 



GM \ ~ 



2c 2 1 



(13) 



where M is the mass and G is Newton's con- 
stant. In these coordinates, the black hole event 
horizon is located at |x| = \ ( ^r ■ Note that, 
because the refractive index becomes infinite 
there, the event horizon is at infinite optical 
distance. A closer examination reveals that as 
one approaches the event horizon, the optical 
geometry approximates more and more closely 
that of hyperbolic space near its boundary at 
infinity as described above. We have recently 
pointed out that this is a universal feature of all 
static ( "non-extreme" ) event horizons and used 
it as a quantitative tool for discussing some of 
their puzzling physics [7 . 

Ricmann also invented an even more gen- 
eral form of geometry, taken up by Paul Finsler 
(1894-1970) called Finsler Geometry in which a 
more general expression replaces (11) and this 



also arises in the optics and acoustics of mov- 
ing media as we shall discuss presently. In the 
meantime we wish to say more about hyper- 
bolic geometry, and in particular the hyperbolic 
plane H 2 . 
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Figure 6: A tiling of the Poincare disk. The Figure 7: The same tiling, this time shown 
edges of the triangles are geodesies. on the upper half-plane model of hyperbolic 

space. 

Source: Three Dimensional Geometry and Topology, Thurston 



2.2 The Hyperbolic Plane 

The geometry of the hyperbolic plane is in- 
timately connected with that of the complex 
numbers. Their introduction dramatically sim- 
plifies many formulae. To see this we adopt 
units in which the radius of curvature UqR is 
set to unity. We further set x = 2R(xi, X2, £3), 
no = 1 in ([9]). The Hyperbolic Plane , H 2 , is 
obtained by setting X3 = . In order to exploit 
the complex numbers we define z = x\ + 1x2 to 
get 



ds 2 



dxl 



dx\ 



\dz\ 



(\-x\-x\Y (1-M 2 ) 2 



(14) 



In this representation of the hyperbolic plane, 
the Poincare Disc, H 2 occupies the interior of 
the unit disc \z\ = 1 and one may verify that 
its geodesies are circular arcs which cut the unit 
circle at right angles. A tiling of H 2 by triangles 
whose edges are geodesies is shown as Figure [6j 
The triangles are all similar, so one can see how 
the apparent length of a measuring rod shrinks 
as we approach the boundary. 

We are free, in addition, to perform a co- 
ordinate transformation to an equivalent form. 
We choose a complex coordinate w = u + iv re- 



lated to z by a fractional linear transformation 

1 + z . . 

w = i- — - . (15) 

This maps the unit disc into the Poincare Up- 
per Half Plane, v > 0. The centre of the unit 
disc maps to the point w = i and the boundary 
of the unit disc to the real axis v = 0. In these 
coordinates the line clement becomes 



ds 2 



\dw\ 2 



du 2 + dv 2 



(16) 



Because, as is easily verified, fractional linear 
transformations take circular arcs to circular 
arcs, the geodesies are again circular arcs. Be- 
cause holomorphic maps are angle preserving, 
the geodeiscs are in fact semi-circles orthogonal 
to the real axis. Figure [7] shows the upper half 
plane tiled with the same geodesic triangles as 
in Figure [6j 

Let's now think of u as horizontal distance 
and v as vertical height in an horizontally strat- 
ified medium in which the refractive index n 
decreases with height. If we assume that over 
a certain range of heights, v, to a good approx- 
imation 

ncx- (17) 
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Figure 8: The hyperbolic plane as the pseudo-sphere in Minkoswki space. Also shown are the 
light-cone X • X = and the plane X° — 0. 



the rays will be semi circles. In this way we 
can easily explain the mirage mentioned in the 
introduction, that in arctic or antarctic regions 
light is bent over icebergs and ships behind ice- 
bergs are observed to float upside in the air 
above. To account for the mirages seen in the 
desert or on hot days driving on motorways in 
which trees or cars seem to be reflected in pools 
of still water, it suffices to take the complex 
conjugate and work in the lower half plane. 

Of course many laws of horizontal stratifi- 
cation n = n(v) will give qualitatively similar 
results, but there is a considerable economy to 
be made by adopting the law (17). If we do 
then the line element ( 16 ) is invariant under 
all fractional linear transformations taking the 
upper half plane into itself. These are of the 
form 



aw 



ad — be = 1 , 



cw 



(18) 



where a, b, c, d are real. This defines the three 
dimensional group SO(2, 1) = SL(2,M)/Z 2 iso- 
morphic with the Lorentz Group of three di- 



mensional Minkowski spacetime. This is no 
accident. Much as it is convenient to think 
of the usual 2-sphere as the set of points at 
a fixed distance from the origin in Euclidean 
3-space, there is a similar interpretation of hy- 
perbolic space as a pseudo-sphere in Minkowski 
space. We consider the space K 3 , but instead 
of the usual dot product, we endow it with the 
Minkowski product: 

X- Y = -X°Y° +X 1 Y 1 +X 2 Y 2 (19) 

where X = (X^^X 1 ^ 2 ) and similarly for 
Y. For a vector with X • X > 0, which we 
call spacelike, we can define the length to be 
|X| = \/X • X. The hyperbolic plane of radius 
R is the set of points defined by 

X-X = -i? 2 , X° > 0. (20) 

The pseudo-sphere is sometimes referred to as 
the mass-shell. This is because when we inter- 
pret the Minkowski spacetime as the geometry 
of special relativity in two spatial dimension^ 
the vectors with P • P < represent momentum 
vectors of particles whose rest mass, m, is given 



In this case we work in units where the speed of light is 1 and a particle of energy E and momentum (pi,P2) 
in the two spatial dimensions would have momentum vector P = (E,pi,p2). 
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byP P 



The mass-shell thus represents 



all the possible velocities for a particle of a given 
rest mass. It is possible to check that any vec- 
tor tangent to this surface is spacelike, so that 
the Minkowski inner product allows us to define 
the length of such a vector. The geometry of 
this surface with this definition of length is that 
of the hyperbolic plane. A Lorentz transfor- 
mation leaves both the Minkowski metric and 



the condition ( 20 ) unchanged and so represents 
an isometry of the hyperbolic plane. To re- 
cover the Poincare disk model, we stereograph- 
ically project the pseudosphere from the point 
( — 1,0,0) onto the plane X° = as shown in 
Figure [8j 



3 Zermelo's Problem and 
Finsler's geometry 

3.1 Zermelo's navigation problem 



In the previous section, we considered the prob- 
lem of finding the ray paths of a sound wave 
propagating in a static medium. We showed 
how Fermat's principle of stationary time leads 
us to consider the geodesies of a Rieman- 
nian metric. We would now like to consider 
how sound waves propagate through a moving 
medium. Fermat's principle continues to apply, 
however we must work a little harder in order 
to translate the principle into a mathematical 
statement. 

We will start off by considering a problem 
proposed by Ernst Zermelo in 1931. Suppose 
a boat, which can sail at a constant rate rel- 
ative to the water on which it sits, wishes to 
navigate from point A to point B as quickly 
as possible. If the water is at rest, then the 
captain should steer along a straight line join- 
ing A to B. More generally, the captain should 
steer along a geodesic if the surface of the wa- 
ter may not be taken to be flat, for example if 
the points A and B are far enough apart that 
the Earth's curvature should be taken into ac- 
count. The navigation problem for the captain 
in this case corresponds to finding the geodesies 
of some Riemannian metric. 



Now let us suppose that the body of water 
on which the boat sits is not at rest, but instead 
moves with some velocity W(x), which we call 
the drift. The absolute velocity of the boat is 
now the sum of two components: the motion of 
the boat relative to the water, v, and W. In 
order to find the fastest route between A and B 
the captain must clearly take the drift into ac- 
count. In order to do this, let us first consider 
the simplest situation, where the surface of the 
water is a plane and the drift W is a constant 
vector, not changing from point to point. This 
situation is shown in Figure [9j We will assume 
for convenience that the speed of the boat rel- 
ative to the water is c, a constant. We will also 
assume that the speed of the drift is less than 
the speed of the boat relative to the water, i.e. 



W 2 < 



Let's define d = AB and work out 



how long it will take the boat to get from A 
to B assuming that v is constant. In this case, 
the position of the boat relative to A at time t 
will be 

x = <(v + W). (21) 

Supposing that the boat arrives at B at time 
T, we have 

d = T(v + W). (22) 

We wish to solve for T as a function of d and 
W. We can do this easily by looking at the 



equations we get from dotting both sides of ( 22 ) 
with W and d: 

d-W = T(v-W + W 2 ), (23) 
d 2 = T 2 (v 2 + 2v- W + W 2 ). 

Noting that v 2 = c 2 and eliminating v • W be- 
tween the equations gives 



T[d] 



(d-W) 



d-W 



C 2_ W 2 ( C 2_ W 2)2 C 2_ W 2- 

(24) 

The condition that the speed of the drift is less 
than c guarantees that the denominators are 
positive and that T[d] > 0, with equality only 
when d vanishes. We notice that if A > 0, then 
T[Xd] = AT[d]. We can also check that T obeys 
a triangle inequality: 



T[d 1 +d 2 ]<T[d 1 ]+T[d 2 }. 



(25) 



This means that the time T which we have 
found is in fact the least time to travel between 
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Figure 9: Zermelo's navigation problem with a uniform drift 



A and B, because we cannot reduce the time 
by travelling along the sides of a triangle with 
base AB. A simple limiting argument shows 
that any curve from A to B will take longer to 
traverse than T. 

Now we can consider a more general prob- 
lem, where the boat is navigating in a Rieman- 
nian manifold M, with metric h — hijdx' '"dx 3 . 
The drift is a vector field W l (x) on M, which 
may vary from point to point and whose length 
is always less than c. Suppose the captain steers 
the boat so as to travel along a curve 7(s), with 
7(a) = A, 7(6) = B. In order to find the time 
taken to traverse this curve, we can approxi- 
mate it with lots of straight sections on which 
the metric and the drift are roughly constant. 
We have worked out how long it takes to travel 
along such a line segment and we can simply 
add these times up. Passing to a limit, we find 
that the time taken to travel from A to B along 
this curve is 



T[l] = \ I F[-'(/M '(.,>/, ( 2 (i, 



where 



F[x\y l ]/c = ^a{x\t)+P{x\f) 
!>,,!/■ !/ J (/'„..'/' 11 - 

C 2_ W 2 + ( c 2 _ ^2)2 



a 



c 2 - W 2 ' 



(27) 



with W 2 = hijW % Wi This problem seems 
somewhat artificial as imagining a boat navi- 
gating a general curved space is rather strange. 



This set up is very natural, however, in the con- 
text of sound rays. If hij is the acoustical met- 
ric of a material (so that a ray moving in the 
direction of the unit vector n z moves with ve- 
locity cj \f hijn % ni) and W l is a bulk motion of 
the material, then T tells us the time the sound 
would take to travel along 7. We are now in a 
position to make a mathematical statement of 
Fermat's principle for a sound wave propagat- 
ing through a moving medium. The sound rays 
are the paths which extremise the time along 
the path, or equivalently the optical length, 
L = Tc: 

5L = cST = 0. (28) 

Before we discuss what this means for the study 
of sound waves in a moving atmosphere, we will 
first discuss briefly a larger class of problems 
which have a similar form to ours. 

3.2 Finsler and Randers geome- 
try 

We have thus far been speaking somewhat 
loosely about geometries, without describing 
exactly what we mean. For our purposes, the 
particular type of geometry which is of interest 
is a Finsler geometry. Although named after 
Paul Finsler, the concept of a Finsler geometry 
was introduced by Riemann in the same lec- 
ture that he proposed what is now known as 
Riemannian geometry. The defining feature of 
a Finsler geometry is that for suitably well be- 
haved curves, one can define a curve length. In 
order to do so, one first has to define a Finsler 



function. The function F of (27) is a special 
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case. A Finsler function, in addition to an as- 
sumption on its smoothness, is required to have 
three properties 

1. Positivity: F[x, y] > 0, with equality only 
ifj/ = 

2. Homogeneity: F[x, Ay] = XF[x,y], for 
A > 

3. Subadditivity: F[x,y\ +2/2] < ^[£,2/1] + 

Given a Finsler function f, we can then dchnc 
the length of a curve from 7(a) to 7(6) by: 

L[l]= I F[f(s),f( S )}ds. (29) 

J a 

Condition 1 ensures that the length of any non- 
trivial curve is positive. Condition 2 ensures 
that the length of a curve does not depend on 
its parameterisation. Condition 3 is necessary 
so that the problem of finding curves of minimal 
length is well posed. This means that we can 
talk about the geodesies of F as being curves of 
minimal length between two points. Note that 
the 'length' defined in this way by the F of the 
previous section is not the usual length of the 
curve, so the geometry defined by this F differs 
from the Euclidean geometry of the plane. 

Whilst Riemannian geometry is fairly well 
understood, Finsler geometry in general is 
much less well studied. This is mainly because 
of the sheer variety of possible Finsler func- 
tions, which can be very exotic. One of the 
simplest classes of Finsler function, into which 
our function F of the previous section falls, are 
the Randers metrics. The Finsler function of a 
Randers metric is given in terms of a Rieman- 
nian metric an and a one-form b; as 



F 



V] = \J a t]{x)y l y3 + h(x)y l 



(30) 



This is a good Finsler function provided 
a l:, bibj < 1, where a 13 is the matrix inverse of 
dij. One reason to study Randers metrics be- 
comes apparent when we consider the equations 
satisfied by a curve which is a critical point of 
(29) when ay is the flat metric in R 3 . Since we 



assume that for the curve x(s) we have x 2 = 1 
and we then find that x(s) is a critical point 
when 

x = xx(Vxb). (31) 

We see that x(s) follows the path of a particle 
of mass m and charge e moving in a magnetic 
field B = m/e(V x b) with unit speed. For 
this reason, the extreme curves of L are some- 
times referred to as magnetic geodesies. For a 
general a^ , we find the generalisation of the 
Lorentz force law in a curved space, with bi 
acting as the vector potential for the magnetic 
field. By extending the notion of a metric to 
allow Finsler geometries, we have brought the 
problem of charged particles in a magnetic field 
into the realm of pure geometry. 

3.3 A uniform magnetic field in 
the Hyperbolic Plane 

As an interesting example, we take a: 1 = u, 
x 2 = v, with v > 0, and we will consider the 



Randers metric given by ( 30 ) with: 

^ 2 a 

61--, b 2 =0. 



- fix 



(32) 



Where p > and a are fixed real numbers. 
This is a good Finsler metric, provided that 
\ct\l p < 1. If a — 0, we see that this Randers 
metric in fact corresponds to the Riemannian 
metric considered in |2.2| that of the hyper- 
bolic plane of radius p in its 'upper half-plane' 
form. If a is non-zero, b corresponds to a vector 
potential for the magnetic field which is every- 
where directed straight out of the plane and 
which has strength B = a/p, independent of 
position. Thus the geodesies of this Randers 
metric will correspond to a charged particle 
moving in a uniform magnetic field in the hy- 
perbolic plane. In order to find the geodesies, 
we have to find curves (u(s),v(s)) which ex- 
tremise the length 



L = p [ ds NE±^ + ™\ (33) 

J \ v pv J 

Such a curve must satisfy the Euler- Lagrange 
equations: 



are free to choose the parameterisation, we can 



= 



ds XvVv 2 + u 2 



\Jv 2 + u 2 a 



p v A 
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Figure 10: Some geodesies of the Randers metric given in Q, with a/p = l/y/2. 



= 



u a 1 

ds \ V y/v 2 + u 2 p v 



(34) 



In the case where a = 0, we know from above 
that the geodesies are circles which meet the 
line v = at right angles. We also know that 
when we have a constant magnetic field in the 
usual Euclidean plane that the particle trajec- 
tories are circles. It seems reasonable then to 
guess that with a non-zero, the geodesies re- 
main circular. We can consider a possible solu- 
tion of the form 



u(s) = uq + r cos s 
v(s) — Vq =F r sins. 



the v = axis are also geodesies, provided 
again that they are traversed in the correct di- 
rection. A little more work shows that in fact 
any geodesic of this Randers metric is one of 
these curves. Figure [10] shows examples of the 
various cases for B = 1/V2- For B < 0, we can 
simply reverse the sense of the curves. 

4 Sound in a wind 



We noticed in §3.1| that Fermat's principle tells 
us that a sound ray in a medium with acoustic 
metric hij with a wind W l will move along a 
(35) geodesic of a related Randers metric defined by 



Notice that the circle is traversed in a clock- 
wise or anti- clockwise direction respectively for 
the two choices of sign. For a general Finsler 
metric, unlike for a Riemannian metric, the di- 
rection of travel is important and a curve will 
only be a geodesic when traversed in a partic- 
ular direction. Substituting into (34) we find 



that we can satisfy the equations provided 



vn a „ 
— = - = ±B. 



(36) 



Thus for < B < 1 we can either have clock- 
wise circles with centres above v = or anti- 
clockwise circles with centres below v = 0. 
Since B < 1, we can interpret both cases ge- 
ometrically as meaning that circles which meet 
the line v = at an angle 8 = cos -1 B are 
geodesies, provided they are traversed in the 
appropriate sense. Taking a limit where r, vq 
get larger and larger with their ratio fixed, we 
find that straight lines making an angle 9 with 



a 



''J 



WiWj 



c 2 -W 2 (c 2 - W 2 ) 2 : 



w 2 ' 



(37) 



where W 2 = h %j W l W^ and W, = h l3 W j . 

Let's firstly consider what this means in the 
case where the speed of sound is constant and 
equal everywhere to c and the wind speed W 
is small in magnitude compared to c. This is 
a reasonable approximation for sound waves in 
a realistic atmosphere. Typically W/c < 0.03 
and even in the strongest hurricanes, W/ c does 
not exceed 0.3. For a constant, isotropic, speed 
of sound, the acoustic metric is simply 



(38) 



If we work to first order in w/c, then we find 
that 

1 W 

"■j h = ~^a- ( 39 ) 
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Making use of (31 1 and keeping track of the fac- 



tors of c, we find that the path followed by a 
sound ray is that of a particle of mass m and 
charge e moving with speed c in a magnetic 
field given by 



B 



-(V x W) 



(40) 



where lu is the vorticity of the wind. This jus- 
tifies our assertion in equation ([I]) that the vor- 
ticity of the wind acts like a magnetic field on 
the sound rays. 

A simple consequence of this correspon- 
dence is observed by seismologists measuring 
the oscillations of the Earth after a large earth- 
quake. The spectral peaks are split by the 
Earth's rotation. From our magnetic point of 
view, the effect of the rotation is to give rise 
to a constant magnetic field inside the Earth. 
The splitting of the spectrum is in precise anal- 
ogy with the Zeeman effect which gives rise to 
a splitting of spectral lines for atoms in a mag- 
netic field. 

4.1 A stratified example 

An interesting problem to study involves com- 
bining a varying speed of sound with a wind 
in a stratified atmosphere. We can gain some 
insight by investigating a particular choice of 
acoustic metric and wind. Let's suppose that 
the acoustic metric for the medium at rest is 
the hyperbolic metric discussed previously 



ds = hijdx l dx 3 



, dx 2 + dz 2 



(41) 



We'll take this to model sound rays in an at- 
mosphere varying with temperature near the 
ground, which we assume to be at z — I. The 
speed of sound at ground level is given by c. 
We'll suppose there's a horizontal wind with 
strength proportional to height, so that 



W 



(42) 



Where w is some fixed parameter with units of 
velocity. We can easily calculate the associated 
Randers metric and find the time to go along a 



curve (x(s), z(s)). It is best expressed in terms 
of new variables: 



iVd 2 ' 



(43) 



This just represents a rescaling of the axes both 
parallel and perpendicular to the ground. In 
terms of u, v, the time is given by 



T = 



I 



Vu 2 +' 



(44) 

Thus, the sound rays follow the geodesies of a 
Randers metric of precisely the form we consid- 
ered in £3.3 with B given by w/cl 



As a simple application, we can consider the 
problem of tracing the sound emanating from 
a point source, P at ground level, which is at 
v = 1. We know that the sound rays in the u, v 
coordinates are circles (and lines) which meet 
the line v = at an angle cos" 1 (w/c). We 
also know that these circles have a clockwise 
sense if their centre is above v = and an anti- 
clockwise sense if it's below. This is already 
enoug h in formation to draw the rays shown in 
Figure 11 for the particular value w/c= l/v2. 
We sketch outward directed geodesies which 
have a common starting point, P. The straight 
line geodesic through P plays an important role 
as a separatrix, which separates two different 
behaviours for the trajectories. In this case, 
any ray emitted to the left of the separatrix 
will return to the ground to the left of P, while 
a ray emitted to the right returns to ground 
to the right. If we suppose that P is emitting 
sound uniformly in all directions, we see that 
more of the energy emitted is absorbed by the 
ground on the left of P than on the right. In 
this case three times as much, which we can see 
by considering the angle the straight geodesic 
makes with the ground. 

This example has three free parameters: c, 
w and I. By choosing these appropriately, we 
can match (for example) the speed of sound, 
the speed of the wind and the the wind shear 
at ground level to a real, more complicated pro- 
file. In fact, in jS] we showed that it is possible 
to construct a model based on the hyperbolic 
plane with four free parameters, so that one can 
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Figure 11: Some sound rays emanating from a point on the ground for w/c = 1/^/2,1 = 1. 
The wind velocity at various heights (including the hypothetical extension below ground level) 
is shown at the left of the diagram. The rays are shown extending as dashed lines below ground 
level to show that they meet z — at an angle of 7r/4. 



additionally set the rate of change of speed of 
sound with height at ground level as well. 

5 Conclusion 

In this article we have sought to show that the 
motion of a charged particle moving in two- 
dimensional Lobachevsky space, or the hyper- 
bolic plane, equipped with a non-uniform mag- 
netic field can provide a useful model for sound 
rays in a moving medium with a gradient in 
the refractive index. We have mentioned that 
in its three-dimensional version the geodesies of 
Lobachevsky space provide a useful model for 
the motion of light rays near the event horizon 
of a non-rotating black hole. If the black hole 
is rotating then Coriolis type effects, referred 
to in General Relativity as the rotation of iner- 
tial frames, provide an effective magnetic field. 
These two examples by no means exhaust the 
possible applications of hyperbolic geometry to 
physics. 

Two-dimensional surfaces abound in nature 
and if they have negative Gauss curvature, 
sometimes called anti-clastic at a point, the sur- 
face cannot lie on one side of its tangent plane 
at that point. Thus a finite smoothly embedded 
surface without edges in Euclidean space can- 
not have everywhere negative Gauss curvature, 
but a finite portion of a surface with edges may. 
A simple example is provided by a holly leaf. 
An example of great current physical interest, 



following the 2010 Nobel Prize to Andre Geim 
and Konstantin Novoselov is a graphene surface 
containing topological defects called disclina- 
tions in which some of the hexagonal lattice 
cells have been replaced by heptagons. The 
electrical and other properties of such surfaces 
are of great interest, and their study entails 
solving the Dirac equation in a portion of two- 
dimensional Lobachevsky space. The motion 
of charged particles on abstract finite Riemann 
surfaces with no boundary or edges which have 
constant negative curvature and uniform mag- 
netic field are of interest in statistical mechan- 
ics since for week magnetic field the motion is 
chaotic or ergodic as it is known technically. 
However as the magnetic field strength is in- 
creased there is a sudden phase transition and 
this ceases to be the case. 

Three dimensional Lobachevsky space has 
been invoked to model some aspects of quan- 
tum dots and the physics of four and five di- 
mensional Lobachevsky space and their confor- 
mal boundaries are currently of intense interest 
by String Theorists since Juan Maldacena sug- 
gested the famous AdS/CFT correspondence 
which has led to a number of break throughs 
in quantum gravity and the quantum theory of 
black holes. Without going into technical de- 
tails, it may be of interest to outline some fea- 
tures of this fascinating idea. Both in Quantum 
Field Theory and in String Theory it is custom- 
ary to work in imaginary time . Thus if we start 
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in Minkowski spacetime with spacetime metric 

ds 2 = -c 2 dt 2 + dx 2 , (45) 

we can pass to Euclidean space with positive 
definite metric 

ds 2 = c 2 dr 2 + dx 2 . (46) 

by setting 

t = it , r real . (47) 

Often calculations may be performed more eas- 
ily in Euclidean space. We then pass back to 
Minkowski spacetime by setting 

t = -it, t real . (48) 

This process is called a Wick Rotation and it 
also works for some curved spacetimes. A case 
in point is Anti-de-Sitter spacetime. This is 
a solution of Einstein's equations with a neg- 
ative cosmological constant. It may be ob- 
tained from Lobachevsky space by a simple 
Wick Rotation. Maldacena's brilliant conjec- 
ture is that there is a precise correspondence 
between String Theory in Anti-de-Sitter space- 
time on the one hand, and a special type of 
quantum field theory called a Conformal Quan- 
tum Field Theory, on the other hand, the lat- 
ter being defined on the conformal boundary 
of Anti-de-Sitter spacetime. The conformal 
boundary of Anti-de-Sitter spactime is confor- 
mally related to Minkowski spacetime If we 
"Wick rotate" this conjecture we are led to con- 
jecture a correspondence between String The- 
ory in Lobachevsky space and Quantum Field 
theory on its conformal boundary, the latter 
being conformally related to Euclidean space. 

We hope that in this article we have made it 
clear that not only is a knowledge of hyperbolic 
geometry and Lobachevsky space useful for un- 
derstanding trafffic noise, but it has a much 
much wider range of applications in theoreti- 
cal physics; from cosmology to condensed mat- 
ter physics to String Theory and Planck scale 
physics. We commend to the interested reader 
its study and further exploitation. 
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